HEPHY-PUB 733/00 
UWThPh-2000-3 8 
bcp-ph/0010078 
September 2000 



INSTANTANEOUS BETHE-SALPETER EQUATION: 
(SEMI-) ANALYTICAL SOLUTION 

Wolfgang LUCHAQ 

Institut fur Hochenergiephysik, 
Osterreichische Akademie der Wissenschaften, 
Nikolsdorfergasse 18, A-1050 Wicn, Austria 

Khin MAUNG MAUNGQ 

Department of Physics, Hampton University, 
Hampton, VA 23668 

Franz F. SCHOBERLfl 

Institut fur Theoretische Physik, Universitat Wien, 
Boltzmanngasse 5, A-1090 Wien, Austria 

Presented by F. Schoberl at the International Conference "Quark Confinement and 
the Hadron Spectrum IV, " July 3-8, 2000, Vienna, Austria 



Abstract 

The Bethe— Salpeter equation for bound states of a fermion-antifermion pair in the 
instantaneous approximation for the involved interaction kernel is converted into an 
equivalent matrix eigenvalue problem with explicitly (algebraically) given matrices. 
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1 The Instantaneous Bethe Salpeter Equation (IBSE) 

For a system of massless fermion and antifermion forming bound states with the 
"pion-like" spin, parity, and charge conjugation quantum numbers J PC = 
the (homogeneous) Bethe-Salpeter equation, in free-propagator approximation 
and instantaneous approximation for the involved interaction kernel, reads fcta 
time-component Lorentz vector interaction (i.e., the Dirac structure 7°£g)7°)IIra 

oo 

-^-^V (k,k')^ 2 (k')=M^ 1 (k) , 

o 

oo 

/c\k' k' 2 
—^ r V 1 {k,k')* 1 {k')=M* 2 {k). (1) 

o 

In this set of coupled equations for the two relevant radial Salpeter amplitudes 
and ^2 in momentum space, with the bound-state masses M as eigenvalues, 
the interaction potential V(r), usually formulated in configuration space, enters 
in form of its standard Fourier-Bessel transforms Vi,(fc, k'), L = 0,1. We adopt 
a linear potential V(r) — Xr (A > 0) as a simple model for quark confinement. 
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2 Efficient Method of Solution: Expansion in Terms of Basis States 

By insertion of the first of Eqs. (|l|) into the second and by expansion in terms of 
sets (distinguished by the angular momenta i = 0, 1) of basis states for L 2 (R + ) 
— with configuration and momentum-space representations (jyf (r) and <j>\ (p 
resp. — the solution of the IBSE ([!]) reduces to the diagonalization of the matri: 

oo oo oo 

M l3 = 4 f dk fc 4 $°> (fc) 4>f (fc) + 2 f dk k 3 $ 0) (fe) / V (k, k') 4>f (fc') 



dfefc 2 $ 0) (fc) I ^- V1 (k,k')^\k') 







dfcfc 2 $ 0) (fc) I ^- Vl (k,k') f^lv Q (k',k")^\k") . (2) 



(2tt) 2 1V ' 'J (2tt 



Allowing these basis functions to depend on a variational parameter /i > gives 
us more freedom in the search for solutions of the IBSE. All integrations in M.m 
are evaluated by (truncated) expansions, with the (/i-independent) coefficientsB 

oo 

^ = ■^1 Mk 4 4> { ?{k)4>f(k) , (U.2.... , 



oo N 

ba = - [ dk k 3 $ 0) (fc) 0f (fc) , k $ 0) (fc) = /x X) b ^ (*) > 





oo 



iV 



dfc fc 2 #w (k) 0<°> (fc) , ^ (0) (fc) = 2 CJj 0f (fc) , 

OO jy 

da = - [ dk k 3 (fc) (fc) , fc $ 0) (fc) = n J2 d ^ $ x) . 





oo 



3=0 
N 

V$ = vl j drr*4>f\r)<i>f{r) , r ^ £) (r) = - £ K.f <^f (r) , £ = 0,1 



The explicit algebraic expressions of all these matrices may be found in Ref. 3.[] 

"Let's mention a numerical problem noted for Mathematica 4.0: for, e.g., the matrix element 
^49, 49 = 101/(2V 689) = 1.924 Mathematica finds exactly this value for a working precision of 
40 digits but the nonsense value —1.675 X 10 22 for the default working precision of 16 digits. 



In this way, the IBSE ([[]) is converted to an eigenvalue problem for the matrixt 

N N N 

Ma = 4 M 2 l\f + 2 A J2 b ri V r f +2X^,^2^ d v V ™ 



r=0 s=0 



' r=0 s=0 t=0 

3 Analytical Results (for Both Massless and Massive Constituents) 

For a matrix size less than or equal to 4, the diagonalization of the matrix Mij 
may be even performed analytically. In the one-dimensional case, we find, aft 
minimizing w.r.t. the variational parameter for the lowest bound-state mass 



/2A 

Fori A = 0.2 GeV 2 , this expression gives M — 1.696 GeV, only 2.4 % away from 
the numerical result M = 1.656 GeV, obtained for 15x15 matrices and N = 49. 
For a nonvanishing mass m of the bound-state constituents, we get accordingly 

,,2 o 2 8896 , 23 / 128A \ 2 . , n . 

M 2 = 8m 2 + —- A+- m^O . 

315 7r 7 \45 nm J 

4 Relations Between Matrix Elements and Accuracy of Expansions 

Our final question concerns the errors induced by the necessary truncations of 
the expansion series. The expansion coefficients bij , Cy , dij are not independent 
but should satisfy (clearly, only in the limit N — + oo exact) relations of the kind 

N N N N 

^ t c ri c rj = 3ij ! ^ ] c ri ^rj — ^ f d ri C r j — bij , ^ ] d ri d r j — I- . 
r—0 r—0 r—0 r—0 

For 15 x 15 matrices and TV = 49, these relations are fulfilled with relative errors 
less than 3 %. For comparison, some integrals in (|2|) may be evaluated exactly.0 
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